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The  Effect^of  Strain-hardening  on  the 
D e for  mat  i on  _of_  an  Annular  ..Slab^- 
by 

P.  G.  Hodge,  Jr.2 

A procedure  is  outlined  for  obtaining  the  stresses  and 
strains  in  a circular  slab  with  a cutout,  subject  to  uniform 
biaxial  teision,  An  arbitrary  stress-strain  curve  in  tension 
is  approximated  by  any  number  of  straight  line  segments.  For 
biaxial  spates  of  stress  the  material  is  assumed  to  satisfy  a 
flow  law  based  on  the  maximum  shear  stress,  and  to  be  incom- 
pressible throughout.  The  general  equations  are  given  and 
then  simplified  by  assuming  that  boundary  motions  may  be  neg- 
lected i:*  the  strains  are  small,  and  that  elastic  strain  com- 
ponents may  be  neglected  if  the  strains  are  large.  For  the 
case  of  linear  strain-hardening  a complete  solution  is  given 
in  closed  form.  If  the  rate  of  strain-hardening  is  small, 

-t' 

these  results  may  be  further  simplified. 


*The  results  presented  in  this  paper  were  obtained  in  the 
course  of  research  conducted  as  a Consultant  under  Contract 
N7onr-358lO  between  the  Office  of  Naval  Research  and  Brown 
University. 

department  of  Mathematics,  University  of  California, 

Los  Angeles.  Jun  ASMS. 


INTRODUCTION 


Consider  an  annular  slab  of  an  incompressible  material 
of  uniform  thickness  H0,  subjected  to  a uniform  biaxial  tension 
of  magnitude  EH0\on  the  outer  edge  (Fig.  1).  It  is  assumed 
that  the  slab  is  in  a state  of  generalized  plane  stress,  even 
though  the  thickness  may  vary  as  the  deformation  progresses. 

Due  to  symmetry,  the  mechanical  state  of  the  slab  may  be  spec- 
ified in  terms  of  the  following  four  functions  of  position  and 
time  i 


d = <Jr(r, t)/E  reduced  radial  stress, 
s=c^(r,t)/E  reduced  circumferential  stress, 


1 


/( 1) 


u=ur(r,t)  radial  displacement, 

h=h  (r,t)  thickness, 

where  E is  Young's  modulus. 

Regardless  of  the  relations  between  stress  and  strain, 
the  slab  must  satisfy  the  equation  of  equilibrium 

(rhe)  - hs  * 0,  (2) 

and  the  condition  of  incompressibility 


r ^ (rhv)  = 0.  C3) 

Hero  v is  the  particle  velocity: 

V * (1  * I?*  It* 

So  long  as  the  load  is  increasing,  the  material  will  be 
assumed  to  be  elastic  provided  that  the  maximum  shearing 
stress  is  less  than  half  the  initial  yield  stress  in  simple 
tension.  It  may  bo  verified  a posteriori  that  the  principal 
stresses  satisfy 


°«  > «r  Z s o, 

so  that  the  material  is  elastic  provided 


<>0 
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f (d  ,s)  = s-p,<0, 


(5) 


where  p.  is  the  reduced  initial  yield  stress  in  simple  tension. 
If  Inequality  5 is  valid,  then  two  additional  equations  are 
furnished  by  Hooke's  law: 

er  = 5u/9r  - (2dr-dQ)/2E  = d-s/2,  1 


eg  = u/r  = (2d0-dr)/2E  = s-d/2. 


? (6) 


The  plastic  domain  is  characterized  by  the  relations 

s-p>0,  §>0. 

The  general  form  of  flow  law  obtained  by  Prager  and  Fodge 
[l]^  may  be  written  in  terms  of  the  reduced  principal  stress 
components 


and  two  similar  expressions  for  e and  e^M.  Fere  e , i 2"  > 

• 

e^'1  are  the  plastic  strain  rate  components,  d^'  , dg'*  d^'  are 
the  reduced  stress  components,  f is  the  yield  function,  and  F 
is  a function  of  stress  to  be  determined  from  the  stress-strain 
curve.  If  f is  given  as  in  Eq.  5,  the  stress-strain  law  becomes 

er"  = 0,  60"  = F s.  (7) 


Let  the  given  stress-strain  curve  be  approximated  by  a series 
of  straight  line  segments,  (Fig.  2),  and  let  ak  be  the  inclina- 
tion of  the  segment  in  the  strain-reduced  stress  plane. 

In  simple  tension,  the  relation  between  plastic  strain  rate  and 
stress  rate  is  then  given  by 


R*x"  = cot  ctfc  - °x^E  ^k+l  » 


Numbers  in  brackets  refer  to  references  collected  at  the  end 
of  the  paper. 
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so  that  Eq.  7 becomes 

er"  = o,  *Q«  = cot  ak  s,  < s £ 4k+1. 


Finally,  since  the  elastic  components  of  strain 
satisfy  the  differentiated  form  of  Hooke's  law, 
stress-strain  relations  are 

er  = 3u/9r  = 5 -s/2, 


1 


/ 


Sq  = u/r  = (1+cot  ak)s-3/2, 


rate  must 
the  plastic 


Equations  2,  3?  and  6 or  8 furnish  four  equations  for  the 
four  unknowns  <3,  s,  h,  and  u.  At  the  boundaries,  the  stresses 
must  be  in  equilibrium  with  the  applied  tractions,  so  that 

<j(a)  = 0,  (9) 

h(d)d(d)  = H0\.  (9b) 

In  addition,  all  stresses  and  strains  must  be  continuous 
throughout  the  slab,  so  that  if  r = pk  is  the  radius  which 
divides  one  analytic  form  of  solution  from  another, 
c(pk)  = c(pk),  s(pk)  = s(pk), 


(10) 


ELASTIC  SOLUTION 


Since  h and  u are  both  unknown,  Eqs.  2 and  3 are  non-linear, 
and  at  present  there  does  not  exist  a general  solution  to  the 
boundary  value  nroblem  defined  in  the  previous  section.  However, 
as  long  as  the  strains  are  small  compared  to  one,  it  appears 
reasonable  to  neglect  motion  of  the  boundaries  and  changes  of 
the  thickness  h.  For  a perfectly  plastic  material [ 2]  it  is 
known  that  the  errors  introduced  by  these  approximations  are  of 
order  p~0.001. 
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Regarding  h as  constant,  3q.  3 can  no  longer  be  satisfied 
exactly,  and  Eqs,  2 and  6 become  linear.  The  solution  for  an 
elastic  material  is  readily  found  to  be 

u = Br  + C/r,  I 


(11) 


2 (B  + 0/3^) 

where  B and  C are  to  be  determined  from  the  boundary  conditions. 

For  X sufficiently  small,  the  slab  will  be  elastic  throughout 

and  the  boundary  conditions  9 may  both  be  applied  to  determine 

the  constants.  Since  boundary  motions  are  neglected,  the 

initial  values  of  a and  d mav  be  used.  The  resulting  solution  is 

2/„2 


u/r  j \ l±3a0  /r 
? = — — * » 


du/dr  f 


2 l-a0^/d0^ 


oL  - s l+a02/r2 
•/ ' ' 


(12) 


The  above  solution  will  be  valid  so  long  as  Inequality  5 is 
satisfied  throughout  the  slab.  Since  s is  a maximum  at  the 
inner  boundary,  the  maximum  load  for  which  this  is  true  is 
found  by  setting  s(a)  * p.  and  solving  for  X : 


X.  k ji  (l-a02/d02). 


(13) 


For  values  of  X somewhat  larger  than  X -p  a plastic  ring  will 
form  against  the  inner  boundary.  If  p denotes  the  radius  be- 
tween the  elastic  and  plastic  regions,  then  the  slab  is  still 
elastic  forp<r<dQ.  In  this  region,  Eqs.  11  are  still  valid, 
but  the  boundary  condition  9a  is  no  longer  applicable.  In 
its  stead,  the  condition 
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s ( p ) = 2(3+C/3p*)  * p 

must  bo  satisfied.  Solving  this  for  G and  substituting  into 
Eqs.  11  one  obtains  the  elastic  solution  in  the  form 


y4  . Bn ; 3 4i  * a..4> 

9u/9r  J + o Jaw  - 


CJ 

s 


» 2 2 
= 2B(1  ± -C)  + p 

r c- 


(1*0 


The  boundary  condition  at  r = d will  be  used  to  determine  a 
relationship  between  \ and  p , once  the  constant  3 has  been 
determined  from  the  plastic  solution. 

ELASTIC-PLASTIC  SOLUTION-FIRST  STAGE 


For  values  of  the  applied  load  slightly  greater  than  X.-^, 
the  plastic  portion  of  the  slab  will  correspond  entirely  to 
the  first  plastic  segment  AB  of  the  stress-strain  curve  (Fig. 
1).  So  long  as  part  of  the  slab  is  still  elastic,  the  plastic 
stress-strain  laws  (Eqs.  8)  may  be  integrated  to  yield 

e _ * au/3r  = o-s/2 

r (15) 

eg  = u/r  = (1+cot  a)s-d/2-p  cot  a. 

The  constants  of  integration  were  determined  from  the  fact  that 
when  a given  particle  first  enters  the  plastic  domain,  the 
strains  must  be  given  by  Eqs.  6 and  s must  equal  p. 

If  Eqs.  15  are  solved  for  s and<j  in  terms  of  u > and  the 
results  substituted  into  the  equilibrium  equation,  the  result- 
ing linear  equation  for  u is  readily  solved  to  yield 


u/r  = Kr11"1  + Lr*’n“1  + k P, 

(16a) 

8u/ar  = nTr11"1  - nLr”*1”1  + ± p, 

(16b) 

„ „ |0E  rn-l  + r-n-l  + n, 

(16c) 
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s 


_ 2n^K  _n-l 
- ‘ 2-n  r 


. 2n2L 
~2 +n 


+ Jk. 


(I6d) 


Here  n is  defined  by 


n = (1+cot  a)“1/2,  (17) 

and  functions  of  a have  been  replaced  by  functions  of  n 
throughout. 

The  four  functions  given  in  the  plastic  region  by  Eqs.  16 
and  in  the  elastic  region  by  Eqs.  14  must  be  continuous  at  the 
elastic-plastic  boundary  r = p.  Together  with  Eq.  9a,  this 
provides  five  equations  to  determine  the  three  constants  B,  K, 
and  L.  However,  in  view  of  the  continuity  imposed  in  integrat- 
ing the  plastic  stress-strain  law,  two  of  these  equations  are 
redundant.  After  some  computation  it  is  found  that 


B * Hi-  1/A)|i, 

K * ^ t pl+n- 

l = Sza  | pl- 


us) 


where 

A = (p/a0)1+n  + (p/a0)1_n.  (19) 

The  substitution  of  Eqs.  18  into  Eqs.  16  and  14  th«n  yields  the 
following  complete  solution: 

Plastic  (aQ<r<p): 

JJ  = k all  + 2+n  (£.)l+n  „ 2^n  <fi.)l-nj  (20a) 

r | nA  r nA  r I 


8u 

8r 

d 

s 


(£)1«-|  (!>1'n]* 


(£  )l+n  n (£  )1- 

r A r 


n . 
» 

_ j 


(20b) 

(20c) 

(20d) 
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Elastic  (p<r£  d0): 

au%?}  - *£  - •1-^2  1 (20o,  f) 

-,fr 

' li  ' * J • (20g,h) 

Equations  20  give  the  complete  solution  in  terms  of  the  posi- 
tion of  the  elastic-plastic  boundary.  The  load  necessary  to 
extend  the  plastic  region  to  a given  boundary  p is  obtained 
from  the  boundary  condition,  Eq.  9b s 


\ = d (d0)  = \i\l  - HP 2/d02 

L A 


(20i) 


Since  the  parameter  n occurs  as  an  exponent  in  Eqs.  20,  the 

preceding  solution  is  not  well  adapted  to  computations.  For 

most  real  materials  n is  small  compared  to  one,  so  that  it 

appears  reasonable  to  approximate  the  solution  and  obtain  a 

more  useful  formulation.  This  may  be  done  by  expanding  the 

solution  as  power  series  in  n.  Since  all  the  quantities  are 

even  functions  of  n,  the  following  expressions  are  correct  up 
L 

to  terms  n i 

Plastic  (a0<r<p): 

u i J"  So  . 0 »o  .P 
^ = 4<|1  +7  + 2 - ldg- 

+ 1?  T [(3+2  log  £)(log  f)  -3(1+2  logf)(log^)1j,(2la) 
« * p|l  - 2 -n2  log  £2  [log  £ + (log  £)2-(log  £-)2l};(21b) 

= p.  4 1 - .12  - \ n2  *fi  [(log  t )2  -(log  £jr]l,  (2lc) 

L r r r aQ  J 

- + l««f  } 


du 

Sr 


(2ld) 


Elastic  (p,<r<d0): 

(1+3  i)  [1  - ^ (log  i-)2  3 V 
r2  2 a0  J 

Q> 

\ — 1 

(21e,f ) 

<*1 

i 

?= nil  - ™ 

(liP£)[l-.n.3.  (log  i.  )2  ])j 
r2  2 a0  J 

(21g,h) 

SJ 

L 2P 

Load: 

X = \i< 

!~i-  iq(i  + £i 

y)  [1-  (log  JL  )?  ] k 

(2li) 

1-  2p  d0' 

? 2 a0  J 

CLASTIC-PLASTIC  SOLUTION -SECOND  STAGE 


The  solution  given  by  Eqs.  20  (or  approximately  by  Eqs.  21) 
is  valid  so  long  as  the  following  two  inequalities  are  both 
satisfied: 

P <.d0»  s(a)l  v • (22) 

As  the  load  is  increased  so  that  one  or  the  other  of  these 
inequalities  is  violated,  a variety  of  behavior  may  be  ex- 
hibited depending  U">on  the  dimensions  of  the  slab  and  the 
parameters  in  the  stress-strain  curve. 

Consider  first  the  case  where  these  quantities  are  such  that 
the  second  Inequality  22  is  violated  while  the  first  remains 
valid,  and  lot  p^  be  the  value  of  p for  which  this  occurs.  To 
the  approximations  used  in  Eq.  21d  it  follows  that  p^  is  given  by 

p (1+n2  logp-j/a0)  = v . (23) 

Since  by  hypothesis  p^  is  less  than  d0  the  parameters  of  the 
problem  must  satisfy 

p (1+n2  log  d0/ac)>  v.  (24) 

The  approximate  load  for  which  this  first  occurs  is  obtained 
by  solving  Eq.  23  for  p^  and  substituting  the  result  into  Eq. 
21i. 
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For  values  of  X somewhat  greater  than  X 2 the  solution  in 
the  slab  will  consist  of  three  regions.  If  £ represents  the 
radius  where  s = v,  then  for  aQ  < r < C the  material  is  plastic 
corresponding  to  the  segment  BJ  of  the  stress-strain  curve 
(Fig.  2);  for  £ < r < p the  material  is  plastic  corresponding 
to  AB;  and  for  p < r < dQ,  the  material  is  elastic,  correspond- 
ing to  OA.  For  the  latter  two  regions,  the  solution  is  given 
by  Eqs.  17  and  l4,  respectively,  where  the  boundary  constants 
must  be  determined  anew. 

For  the  innermost  region,  the  stress-strain  law  Eq.  8 may 
again  be  integrated  to  yield 

er  = 9u/  9r  = <3-s/2  , 


> (25) 

zq  = u/r  * (1  + cot  p)s-o/2  - v (cot  {3-cot  a)  - p.  cot  aj 
As  in  the  previous  section,  those  relations  may  be  solved  for 
<3  and  s and  the  results  substituted  into  the  equilibrium  equa- 
tion. The  resulting  equation  is  readily  solved  for  u: 


whore 


2 

K 


Pr®"1  + Qr"®”1  + k 


1-n^  m2 

13  5 


+ 


V, 


m ss  (1  + cot  p)“1//2. 


(26) 

(27) 


The  constants  P,  Q,  K,  L,  and  3 may  now  be  determined  from  the 
continuity  conditions  (Eqs.  10)  at  C and  p . The  boundary 
condition  9a  then  determines  the  relation  between  p and  C , 
while  Eq.  9b  furnishes  the  value  of  the  load.  The  computations 
and  results  become  quite  complicated  and  are  not  reproduced 
here. 
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Consider  now  the  other  possibility  characterized  by  the 
reverse  of  Inequality  24.  In  tMs  case  the  slab  will  become 
fully  plastic  for  a load  X^  obtained  by  sotting  p equal  to  dQ 


in  Eq.  20i  or  21i.  Thus  fro  i Eq.  21i, 

,2 


X'  = LL 
2 * 


i - sa  + 


5°(lo5  ) 

2 U r\ 


do 
a© 


(28) 


For  values  of  \ slightly  greater  than  Eqs.  16  will  be  valid 
throughout  the  slab.  The  constants  K and  L may  be  determined 
directly  in  terms  of  the  load  X by  substituting  Eq.  16c  into 
Eqs.  <■>.  While  the  results  mav  be  obtained  accurately,  it  is 
somewhat  simpler  to  introduce  the  approximate  power  series  in 
n directly.  After  some  computation,  the  complete  solution  is 
found  to  be 


— _ ^ ^2  (ao/dp)4  " 4 

r " r?  r log  (dQ/a0) 


3u 


dr/r  i aQ  d0  , r i 

2 log  (d0/a0)'  ip*1  (1*2  log  “)  + <WK2  log  — -lj 

+ i I*  > 


'd0/r 

log  <do/a0) 


"4  T*  log  IT  + ^“4  ) loS 

dQ  r r ao 


+ 4. 


(2<?) 


Q = d0  (ao/d0)p  + X -p, 

' T — Yog ' TcfoT^T-  +<  “ * 


whore  only  the  loading  terms  in  n^  have  been  retained. 

Equations  29  will  be  valid  until  s(a)  =v  , assuming  that 
the  strains  may  still  bo  considered  small.  On  the  other  hand, 
assuming  for  simplicity  that  the  entire  stress-strain  curve  is 
approximated  by  the  three  segments  in  Fig.  2,  the  solution 
obtained  from  Eqs.  27,  17,  and  14  is  valid  until  p = dQ.  For 


1 


Bll-18 


12 


slightly  greater  loads  in  either  case,  the  slab  will  be  plastic 
throughout,  but  will  be  governed  by  the  portion  3C  of  the  stress- 
strain  curve  for  a0<r<£;  and  by  AB  for£  ,<  r<  ^0.  Therefore, 
the  solution  will  consist  of  Eqs.  27  for  a0>£r<£and  ^Qs* 
for££  r < d0,  with  the  four  constants  P,  Q,  K,  and  L to  be 
determined  from  Eqs.  9 and  10. 

If  a more  accurate  approximation  to  the  stress-strain  curve 
is  desired,  the  solution  may  be  obtained  in  similar  fashion, 
although  the  resulting  computations  would  become  exceedingly 
laborious. 

FULLY  PLASTIC  SOLUTION-LARGE* STRAINS 

The  solutions  obtained  in  the  preceding  sections  are  valid 
only  as  h may  be  considered  a constant  and  the  boundary  dis- 
placements neglected.  For  an  incompressible  material  these 
conditions  will  both  be  satisfied  provided  that  the  strains  are 
small  compared  to  one  throughout  the  slab.  This  will  generally 
be  the  case  so  long  as  any  part  of  the  slab  is  elastic,  at 
least  provided  a0/dQ  is  not  too  small.  However,  once  the  slab 
is  fully  plastic,  relatively  large  strains  may  occur  and  the 
previous  solutions  would  no  longer  be  applicable. 

Since  h can  no  longer  bo  considered  constant,  the  equations 
to  be  solved  are  non-linear,  and  it  is  necessary  to  make  some 
other  simplifying  assumption  to  obtain  a solution.  The  elastic 
part  of  the  strain  remains  proportional  to  the  stress  throughout 
loading  and  lunco  oc.n  never  become  largo.  Therefore,  if  the 
total  strains  are  large,  it  appears  reasonable  to  neglect  entirely 
the  elastic  component  and  consider  the  strains  as  purely  plastic. 
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This  is  equivalent  to  assuming  that  Young's  modulus  is  equal  to 
infinity,  so  that  tho  clastic  portion  cf  the  stress-strain 
curve  in  Fig.  2 is  replaced  by  a vertical  line.  F^r  sim- 
plicity n f exposition  it  will  be  assumed  that  the  plastic 
portion  of  the  stress-strain  curve  consists  cf  a single 
linear  segment  (Fig.  3)  although  the  method  can  be  extended 
to  more  general  curves  as  in  the  previous  sections. 

In  tho  case  of  small  deformations,  it'  was  not  necessary 
to  distinguish  between  the  Lagrangian  (original)  and  Eulcrian 
(current)  coordinates  of  a particle.  However,  for  the  large 
strains  considered  in  the  present  section  such  a distinction 
is  obligatory.  The  stress-strain  law  is  a function  of  the 
stress  and  strain  of  the  particular  particle  in  quest!:  n, 
and  hence  must  be  expressed  in  terms  of  the  original  coor- 
dinate rD  of  a particle  now  at  r . Thus,  since  clastic 
strain  rates  arc  neglected,  Eqs.  8 are  replaced  by 

e = avii'S  r = o,  (30a) 

e6  = 2 = cot  a § = [ (l-n2)/n2  ] S . (30b) 

The  "dot"  in  Eqs.  30  customarily  refers  to  time.  However, 
since  inertia,  ‘effects  arc  not  considered,  time  may  bo  re- 
placed by  any  convenient  monotonic  function  of  time.  In  the 
present  ease,  the  current  radius  of  the  inner  boundary  is 
such  a function.1*' 


Tho  validity  of  this  stop  may  bo  seen  by  noting  that 
the  rate  at  which  load  is  applied  will  not  affect  tho  sclu- 
tion*  Thus,  in  particular,  tho  lead  may  be  applied  so  that 
tho  inner  boundary  moves  outward  at  a uniform  rate.  Finally, 
tho  unit  of  time  may  be  chosen  so  that  this  rate  is  unity. 
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Therefore  Eq,  30a  states  that 
37-  v (r,a)  = 0. 

While  this  equation  may  easily  be  integrated  to  obtain 

the  velocity,  it  is  mere  convenient  to  first  transform  it 

to  the  Lograngian  coordinate  r0  . Since  the  ’‘time'’  a occurs 

only  as  a parameter, 

22L.  = av/dxo  _ 0 

8r  9r/0ro  u’ 

hence 


9v 

W, 


9 

■Sr: 


< -fe  > = °’ 


■o  — o 

where  u(rc,  a)  is  the  displacement  of  the  particle.  The 
integral  of  this  equation  shov/s  that 

u (rG,a)  = F(r0)  + G(a). 

Since  elastic  strains  are  neglected,  the  displacement  of  all 
particles  is  zero  at  the  commencement  of  plastic  flow.  Thus 


u (r0,a0)  » F(r0)  + G(a0)  = 0. 


FCr^)  can  now  be  evaluated,  the  resulting  displacement  being 

u (rc,a)  = G(a)  - G(aQ). 

In  view  of  the  choice  of  a time  scale,  the  velocity  of  the 
particle  initially  at  the  inner  boundary  (and  hence  always 
at  the  current  inner  boundary)  is  unity,  so  that 


Therefore  G(a)  = a,  and  hence 

u (r0,a)  = a - a0  (31a) 

independently  of  ro  . Further  the  current  position  of  the 
particle  initially  at  r0  is 


/ 
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r(r0,a)  = rQ  + u(rQ,r)  = rp  + a-ap. 

Finally,  tho  velocity  at  all  points  is 

V © 1. 

The  substitution  of  Eqs.  31  into  Eq.  30b  yields 


(31b) 

(31c) 


9s 

da 


n 


r -i? 


hcnco 


r0+a“ao 


(32) 


o(r0,a)  = -A-^  log  (r0+a-aQ)  + f(rQ) 

Tho  function  f (r0)  is  to  bo  evaluated  from  the  continuity 

of  s at  a = ac.  However,  since  there  has  been  a change  of 
assumption  between  the  small  strain  and  large  strain  solu- 
tions, this  continuity  cannot  be  exact.  The  best  one  can  do 
is  to  evaluate  s(rp,a)  so  that  it  is  continuous  with  tho 
solution  given  by  Eq.  20d  for  p =d0.  This  loads,  of  course, 
to  a rather  awkward  expression  for  f(r0).  However,  this  ex- 
pression may  be  simplified  by  neglecting  terms  of  order  n*. 
This  leads  finally  to  _ 

s e p + n2  ML 


log  ZSIEZl S + p £2  log  ifi 
ro  “ 


r o+a-ao 

ro  ro  J 

a result  which  could  also  have  been  obtained  directly  from 


Eq.  21d.  Finally,  since  -a,  is  very  small  compared  to  one, 

p 

the  second  torn  in  tho  coefficient  of  n may  bo  neglected  in 
comparisicn  with  tho  first,  once  the  deformation  has  progres- 
sed, Thus,  to  within  the  approximations  considered, 

s = p + n2  log  JL  , (33a) 

ro 

whore  r and  rQ  aro  related  by  Eq.  31b.  Equation  31d  predicts 
that  s = p throughout  the  slab  at  tho  onset  of  full  plasticity, 
Thus  physically  the  assumptions  arc  equivalent  to  neglecting 
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the  strain-hardening  which  occurs  during  the  elastic-plastic 
stage  of  loading.  In  an  example  considered  later,  this  approx- 
imation appears  to  be  well  justified. 

In  Eq.  31c,  the  functional  forn  of  the  velocity  v is  the 
same  referred  to  either  the  Lagrangian  or  Eulerian  coordinates. 
Thus  the  Eulerian  form  of  Eq.  3 (incompressibility  equation)  may 
be  retained: 


± <rh)  ♦ 
oa 


a 

"5r 


(rh)  = 0, 


where  h is  a function  of  r and  a.  The  solution  of  this  equation 
satisfying  the  condition  h = H0  when  a = aQ  is 

h - <1  - i?-2)  H0  = £°  Hc.  (33b) 

r r 


This  result  could  also  have  been  predicted  directly  as  follows. 

For  velocity  constant,  the  thickness  A r of  a narrow  ring  will 
remain  constant,  hence  by  incompressibility  the  volume  2nr  h (Ar) 
must  also  remain  constant  at  all  times  a. 

The  radial  stress  is  then  determined  by  integrating  the 
equilibrium  equation  (Eq.2)  with  respect  to  the  Eulerian  co- 
ordinate rt 


= JL 

TO 


rrc 


vl  a 


£2  + n2  £2  log  £2  ] dr 
V V T 


[ r-a-(a-aQ)  log  £ ] - £- 
w a ro 

u 


(33c) 


Bll-18 


i »cf 


17 


In  the  above  integrations,  r0  is  defined  as  a function  of  r by 
Eq,  31b.  Finally,  the  load  necessary  to  enlarge  the  inner  radius 
from  a0  to  any  value  a is  obtained  from  Eq.  9b: 

X (a)  » h (d,a)  a (d,a)/HQ 


. u|i  [1  - .5  . (*£2)  log  1] 
d d a 


dr, 


where  the  current  outer  radius  d is  given  by 

d = dg  + a — flg« 

For  a perfectly  plastic  material  [2  ] , it  was  found  that 
the  load  X decreased  as  the  fully  plastic  deformation  proceeded, 
so  that  the  maximum  load  which  the  slab  could  support  was  achieved 
just  as  the  elastic  ring  disappeared.  For  a material  with  infin- 
ite strain-hardening,  on  the  other  hand,  the  load  could  increase 
indefinitely  above  the  fully  plastic  value.  The  question  natural- 
ly arises  as  to  the  rate  of  strain-hardening  which  separates  these 
two  types  of  behavior.  In  other  words,  what  is  the  minimum  strain- 
hardening rate  for  which  the  initial  fully-plastic  configuration 
is  stable? 

Since  X is  given  as  a function  of  a by  Eq.  33d,  this  ques- 
tion is  easily  answered  bv  differentiating  X with  respect  to  a, 
setting  a = aQ,  and  observing  the  sign  of  the  derivative.  If  it 
is  positive,  then  an  increase  in  radius  must  be  accompanied  by 
an  increase  in  load.  After  some  computation,  it  is  seen  that  the 
derivative  at  ao  is  positive  if  and  onlv  if 

n2  > |j,  1 - ap/dp  + log  (d0/a0)  ^ 

log  d0/a0 
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Equation  3^  states  the  condition  that  the  fully  plastic  solu- 
tion be  initally  stable.  As  the  deformation  progresses,  however, 
it  follows  from  Eq.  33b  that  the  thickness  at  the  outer  edge  grows 
less,  so  that  the  same  load  produces  a higher  true  stress.  There- 
fore, for  any  finite  rate  of  strain-hardening,  the  actual  load 
which  the  slab  can  support  will  decrease  after  the  deformation  has 
proceeded  sufficiently  far.  The  value  of  a for  which  this  is  the 
case  can  be  foqnd  by  equating  d\/da  to  zero  and  solving  for  a. 
Although  this  cannot  be  done  in  closed  form,  the  maximum  value  of 
X can  be  found  numerically  for  any  given  values  of  the  parameters. 

EXAMPLE -PERFECTLY  PLASTIC  MATERIAL 

A perfectly  plastic  material  is  characterized  by  a zero  rate 

of  strain-hardening,  i.e.,  by  a * 0.  Therefore,  the  complete 

solution  for  a perfectly  plastic  material  would  consist  of  three 

stages*  fully  elastic,  elastic-plastic,  and  fully  plastic.  The 

elastic  solution  is  still  given  by  Eqs.  12 , while  in  view  of  Eq.  17 

the  other  two  solutions  are  obtained  by  setting  n=o  in  Eqs.  21,  31, 
and  33  respectively*  -\. 

*v 

Elastic-plastic  solution  for  aG  < r < Pi  ~ 1 * 

u/r  * k p [1  + a0/r  + 2 (aQ/r)  log  (p/r)], 

du/9r  = i |i  (1-2  a0/r), 

d = p (l-a0/r),  s = p; 

Elastic-plastic  solution  for  p < r < d0* 

«$*}”  * “ B 1 1 3 f^1’ 

“ [1  • f? T §5 1 * 


r (35) 
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Elr.stic~nlr.stlc  lord: 

X = ,1  [ 1-  - SftSjr  ]. 

Fully  plastic  solution: 

u/r  = (a-a0)/r,  au/ar  = 0,  h = [1-  (a-an)/r  ] hc, 

„ = l r.  n - Cr-an).  log,  (r/r.) 

d s = u . 


I (35) 


X = |t  [l-ao/d0  -(r./do-a0/do)lo2/(do+r.-a0)/a'j>  ]* 

(1  + a/d^  — Sq/^o) . 


In  a previous  paper  [ 2 ] , the  author  obtained  the 
complete  perfectly  plastic  solution  directly  by  a perturba- 
tion method.  It  may  readily  be  verified  that  the  first  term 
of  each  quantity  in  the  perturbation  procedure  agrees  precisely 
with  Eqs.  3r'»^  In  the  discussion  following  the  presentation 
of  [2]  at  the  1952  Annu;  1 Mooting  of  A.S.M.E.  (New  York,  N.Y., 
Doc.  2,  1952),  the  question  was  raised  as  to  whether  or  not 
the  perfectly  plastic  solution  truly  represented  a limiting 
case  of  a strain-hardening  solution;  the  present  example  shows 
that  for  the  problem  considered  it  d^os  indeed. 


The  clastic-plastic  solution  is  given  by  Eqs.  32  of 
[ 2]  , and  the  fully  plastic  by  Eqs.  3°  through  46.  Noto 
that  the  dimensionless  stressos  there  wore  defined  by  divid- 
ing by  twice  the  shear  modulus  2G,  rather  than  Young's  modulus 
£,  and  that  |i  was  defined  by  k/2G.  Since  E = 3G  for  an  in- 
compressible material,  the  results  may  be  easily  converted. 
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EXAMPLE -STRAIN-HARDENING  MATERIAL6 

Some  results  have  boon  worked  cut  for  the  following 
numerical  data: 

p,  = 0.005,  n2  = 0.05,  ac/d0  = 1/3. 

Tho  plastic  part  of  tho  stress-strain  curve  has  been  approx- 
imated by  a single  linear  segment.  Figure  4 shows  tho  load, 
circumferential  strain  u(a0)ac  of  tho  inner  boundary,  and 
reduced  circumferential  stress  s(a0)  at  the  inner  boundary, 
all  as  functions  of  the  radius  "of  tho  clastic-plastic  boundary 
for  the  clastic-plastic  solution  given  by  Eqs.  21.  The  maximum 
strain,  when  the  slab  first  becomes  fully  plastic,  is  seen  to  bo 
greater  than  one  per  cent.  Since  this  strain  can  bo  shown  to 
increase  rapidly  with  further  increase  in  load,  the  fully 
plastic  solution  for  small  strains  given  by  Eqs.  29  docs  not 
appear  to  bo  applicable  for  the  present  example.  Therefore, 
in  Fig.  5 the  load,  maximum  stress,  and  maximum  strain  are 
plotted  as  functions  of  the  radius  of  the  inner  boundary 
according  to  Eqs.  31  end  33  for  the  fully  plastic  slab.-  It 
will  be  observed  that  w’-'ilc  the  slab  becomes  fully  plastic  for 
\ * 0.0034,jt  can  support  a load  X = 0.0085  with  sufficient 
plastic  deformation.  This  is  in  markod  contrast  to  tho  per- 
fectly plastic  slab  which,  in  the  present  example,  becomes 
fully  plastic  for  X = 0.0033  end  can  support  no  greater  lord. 


Tho  author  wishes  to  thank  Mr.  R.K.  Froyd  for  his 
assistance  with  the  computations  of  the  present  section. 
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Figure-  6 is  derived  fr^m  Figs.  4 and  5 and  shows  the 
circumferential  stress  and  strain,  rospoctivoly,  at  the 
inner  boundary  as  functions  of  the  load  for  the  clastic, 
clastic-plastic,  and  fully  plastic  srluticns.  The  discon- 
tinuity between  the  last  two  sections  of  the  curve  in  each 
ease  is  due  to  the  change  of  assumptions  between  the  two 
solutions.  It  is  seen  to  bo  snail  in  comparison  with  the 
other  quantities  involved. 
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